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ABSTRACT 

We present a Chandra study of 38 X-ray luminous clusters of galaxies in the ROSAT 
Brightest Cluster Sample (BCS) that lie at moderate redshifts {z ~ 0.15-0.4). Based 
primarily on power ratios and temperature maps, we find that the majority of clusters 
at moderate redshift generally have smooth, relaxed morphologies with some evi¬ 
dence for mild substructure perhaps indicative of recent minor merger activity. Using 
spatially-resolved spectral analyses, we find that cool cores appear to still be common 
at moderate redshift. At a radius of 50 kpc, we find that at least 55 per cent of the 
clusters in our sample exhibit signs of mild cooling (tcooi < 10 Gyr), while in the 
central bin at least 34 per cent demonstrate signs of strong cooling (fcooi < 2 Gyr). 
These percentages are nearly identical to those found for luminous, low-redshift clus¬ 
ters of galaxies, indicating that there appears to be little evolution in cluster cores 
since z w 0.4 and suggests that heating and cooling mechanisms may already have 
stabilised by this epoch. Comparing the central cooling times to catalogues of central 
Ho emission in BCS clusters, we find a strong correspondence between the detection of 
Ha and central cooling time. We also confirm a strong correlation between the central 
cooling time and cluster power ratios, indicating that crude morphological measures 
can be used as a proxy for more rigorous analysis in the face of limited signal-to-noise 
data. Finally, we find that the central temperatures for our sample typically drop by 
no more than a factor of ~3-4 from the peak cluster temperatures, similar to those of 
many nearby clusters. 

Key words: surveys - galaxies: clusters: general cooling flows - X-rays: galaxies: 
clusters - 


1 INTRODUCTION 

By the time a cluster of galaxies reaches a relatively re¬ 
laxed state, the hot gas in the centre is likely to have a 
radiative cooling time which is shorter than the expected 
cluster age and a temperature which drops toward t he cen¬ 
tre, d ue to energy losses from X-ray emission (e.g.. iFhbiarJ 
Il994f) . This effect, termed a “cooling flow”, provides a mech¬ 
anism by which matter can condense out of the hot intra¬ 
cluster medium (ICM) and is observationally detected as 
an enhanced X-ray surface brightness peak in the cores of 
clusters. The final state of this cooling gas i s still a matter 
of intense debate (e.g., |PetersorLet_alJ| 200l BTamura et alJ 

l200ll : IPeterson et alj[20fl^nKaasto^^STl20o3) ~ as the tem- 

peratures of the cooling gas are found to generally drop by 
less than a factor of ~3-4 with very little cooler gas (pre- 
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sumably due to some form of heating, such as from an ac¬ 
tive galactic nucleus). Despite this cooling “floor”, and the 
fact that the mass of gas cooling to zero represents only 
a tiny fraction of the overall hot phase of the ICM, cool¬ 
ing cores remain an observational phenomenon which inex¬ 
plicably trace the physical processes by which clusters of 
galaxies form. In general, the central regions of a cluster 
of galaxies must remain relatively relaxed both to estab¬ 
lish and maintain a cooling core. Strong disturbances to the 
cluster, such as might be caused by merging for example, 
should mix the cluster gas and could disrupt the cooling 
core. The extent to which the cooling core can be suppressed 
or destroyed, in fact, depends on a number of factors such 
as the initial strength of the cooling core, how off-centre 


McGlvnn & FabiarJ Il984l: lEdee. Stewart. & FabiarJ 

1992 

Allen et alJl20011: IGomez. Token. R.oettieer. & Burna 

2002 

Ritchie & ThomafJ l2002^ . Thus, a proper census of cool- 
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ing cores in clusters as a function of lookback time should 
provide important constraints on the robustness of cooling 
cores, the nature of cluster buildup and allow further refine¬ 
ments to structure formation models. 

Central surface brightness peaks in clusters of galax¬ 
ies (the observational signature of cool cores) appear to 
be remarkably widespread, residing in 70-90 per cent of 
the clusters with 2 < 0.10, and o ccasionally even dominat¬ 
ing their bolometric o utput (e.g., lEdge, Stewart, fc Fabianl 
ll992l:IPeres et al.ll998ti . The ages of these cooling cores, typ¬ 
ically considered to be the time since their last disruption, 
are es timated to be < 4 Gyr in ^50 per ce nt of the ob¬ 
jects iAllen fc Fabianlll997tlAllen et ^l200lll . By observ¬ 
ing galaxy clusters at 2 « 0.15-0.4 (i.e., lookback times 
of 2-4 Gyr), we can probe the conditions of these objects 
during a potentially important disruption period. The ques¬ 
tion of how common cooling cores are in the more distant 
past remains completely unexplored. We report here on a 
large sample of moderate redshift clusters with Chanda ob¬ 
servations. The sub-arcsecond spatial resolution afforded by 
Chandra is extremely important in this effort, as it allows 
the detection and quantification of cooling cores out to much 
larger distances than was previously accessible with past X- 
ray observatories such as Einstein and ROSAT. 

We describe our X-ray sample in analysis of the 
Chandra data in m and finally, our results in m Through¬ 
out this paper, we adopt Hq = 70 km s“^ Mpc“^, SIm = 0.3, 
and Qa = 0.7. Unless explicitly stated otherwise, quoted er¬ 
rors are for a Icr (68 per cent) confidence level. 


2 SAMPLE SELECTION 

We selected ou r sample from the R OSAT Brightest Gluster 
Sample (BGS; lEbeling et al.lll9 9^ an d Extended R OSAT 
Brightest Gluster Sample fEBGS: lEbeling et al.l200nll . both 
of which are derived from ROSAT All-Sky Survey data. 
When combined (hereafter simply BGS), they represent one 
of the largest and most complete X-ray-selected cluster sam¬ 
ples compiled to date. Figure 0 highlights the fraction of 
BGS sources with Chandra observations, shown in the 0.1- 
2.4 keV luminosity versus redshift plane. Importantly, a 
large fraction of the most luminous and most distant BGS 
sources have already been observed by Chandra. We adopt a 
selection criterion of 2 > 0.15 to maximize the average red¬ 
shift of the sample, while still providing a statistically use¬ 
ful number of objects with Chandra observations. Of the 51 
clusters with 2 > 0.15, 38 (75 per cent) have publicly avail¬ 
able Chandra exposures as of 2003 October. The vast ma¬ 
jority of these are of sufficient quality (i.e., >5000 counts) to 
assess crudely the nature of any potential cooling core. Thus, 
our sample should be sufficiently complete to provide reli¬ 
able statistics on luminous clusters of galaxies at moderate 
redshifts. As a local com parison samp le, we use the Bright¬ 
est 55 sample studied bv IPeres et al.l (hereafter B55: Il991i) . 
The B55 sample is a 2-10 keV flux limited sample of X-ray 
emitting clusters which is nearly complete and comprised 
of sources which are all close enough to have been imaged 
properly with previous X-ray instruments (i.e., ROSAT). To 
make the B55 sample more comparable with ours, we have 
cropped two clusters with 2 > 0.15 and excluded clusters 
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Figure 1. X-ray luminosity versus redshift for the BCS sample 
(circles). Filled circles denote BCS clusters observed with Chan¬ 
dra as of 2003 October. Our sample consists of those sources with 
2 ^0.15. 

with Lo.i_ 2.4 keV < 4 X lO'*'*. We have converted the B55 

values to our adopted cosmology. 


3 DATA 

3.1 Preparation 

The Chandra data (see Ta,hle[T} were processed and cleaned 
using the GIAO software and calibration files (GIAO v3.0.2, 
GALDB 2.27) provided by the Chandra X-ray Genter 
(GXG), but also with FTOOLS (v5.3.1) and custom soft¬ 
ware. We began by reprocessing the data with the GIAO tool 
AGISJROGESS-EVENTS to remove the standard pixel 
randomisation and to flag potential AGIS background events 
for data observed in Very Faint (VF) mode. The data were 
then corrected for the radiation damage sustained by the 
GGDs during the hrst few months of Chandra operations us- 
ing the charge t ransfe r inefficiency correction procedure of 
iTownslev et alJ (jioO^).^ Following GTI-correction, we per¬ 
formed standard grade selection, excluding all bad columns, 
bad pixels, cosmic-ray afterglows, and VF-mode background 
events when applicable. We used only data taken during 
times within the GXG-generated good-time intervals. Fur¬ 
thermore, light curves for each dataset were made to screen 
for periods of high background; periods exceeding 3a of the 
mean rate were removed. 

Backgrounds rates were then estimated from local re¬ 
gions as free of cluster emission as possible; in a few instances 
where large clusters were observed with AGIS-S, our back¬ 
ground rate estimates are likely to be contaminated upward 
by faint cluster emission. Ratios of the average 0.3-10 keV 
background count rate for each observation to that of AGIS 
quiescent-background calibration values (i.e., 0.31 cts s“^ 
per AGIS-I chip and 0.86 cts s“^ per AGIS-S3 chip) are 
listed in Table Q The vast majority of the observations have 


^ For details see http://www.astro.psu.edu/users/townsley/cti/ 
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Table 2. General Properties of the Cool Core Sample. 


(1) 

(2) 


(3) 

(4) 

(5) 

(6) 

(■?) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

Cluster 

Position 

(J2000) 


JVh 

kT 

z 

Fx 


Ha 

Mor. 

-Rcent 

^cool 

^50 

*180 

A68 

00** 37"* 06? 1 

-1-09° 09' 

' 34 " 

4.6 

10.0 

0.2546 

5.5 

14.89 

7 

R 

41.3 

24.0+10.6 

24,0+10.6 

57 5+°-° 
^^•^-22.6 

All5a* 

00'*55"*50®9 

-1-26° 24' 

'35" 

5.2 

9.8 

0.1971 

9.0 

14.59 

Y 

D 

17.4 

l' 6 lo .2 


23.0+1“ 

— 4.1 

All5b* 

00'*55"*59?4 

-1-26° 20' 

'13^^ 

5.2 

9.8 

0.1971 

9.0 

14.59 

7 

D 

43.1 

22.0^11 

22 . 0 ;“:^ 

34 1 + ®2.1 
^ -10.4 

A267 

0l'*52"*42?l 

-l-01°00' 

'38'^ 

2.8 

9.7 

0.2300 

6.2 

13.71 

N 

R 

36.9 

^^•-‘^-13.1 

25 1 +^'2 

22 . 0 ;™/ 

A520 

04'*54"*09?8 

-1-02° 55' 

'16^^ 

7.4 

9.8 

0.2030 

8.4 

14.44 

N 

D 

52.7 

14.8+l^g® 

13.6+1“/ 

— 4.4 

24 5 + 23^2 

A586 

07'*32"*20?3 

-l-31°37' 

'56^^ 

5.4 

8.7 

0.1710 

9.1 

11.12 

N 

R 

30.2 

0 q“I"4.3 

Q q+4.3 

15.8;iJ 

A665 

08'*30"*58?0 

-1-65° 50' 

' 54 ^^ 

4.0 

8.3 

0.1818 

11.8 

16.33 

N 

D 

38.4 

18.2+“j;9 

18 . 2 ;// 


A697 

08'*42"*57?7 

-|-36°2l' 

’ 57 ^^ 

3.4 

10.5 

0.2820 

5.0 

16.30 

N 

R 

43.6 


28 . 4 ;“/, 

25 . 3 ;“+“ 

A750 

09'*09"*12?2 

-l-10°58' 

'36'^ 

3.6 

8.1 

0.1630 

8.3 

9.30 

N 

R 

22.9 

n rt + 2.5 
'■^-1.6 

'•^-1.6 


A773 

09'*17"*52?9 

-|-51°43' 

'37^^ 

1.5 

9.4 

0.2170 

6.7 

13.08 

N 

R 

33.1 



63.4j;«4«4 

A781 

09'*20"*25?4 

-1-30° 30' 

'13^^ 

1.9 

10.8 

0.2984 

4.7 

17.22 

7 

2 

81.7 

17.4+«;9 


17.4;i“ 

A963 

10'*17"*03?5 

-1-39° 02' 

’ 54 ^^ 

1.4 

8.6 

0.2060 

5.9 

10.41 

N 

R 

20.8 


4 5+°-® 
^-0.7 

23.1+1“ 

— 4.5 

A1204 

ll'*13"*20?4 

-1-17° 35' 

'38^^ 

1.4 

7.5 

0.1904 

5.0 

7.62 

Y 

R 

9.5 


2.1‘''°-^ 

-0.3 


A1423 

ll'*57"*17?2 

-1-33° 36' 

'39'^ 

1.6 

8.5 

0.2130 

5.3 

10.03 

N 

R 

24.4 


A 9+1-7 
-1.1 

22 . 5 ;// 

A1682 

13'‘06"‘51?8 

-1-46° 33' 

'13'^ 

1.1 

8.9 

0.2260 

5.3 

11.26 

N 

2 

49.4 

34 6+^^ ® 

^ "-18.2 

34 » 

^ ^-18.2 

113 

Al758a 

13'‘32"*43?5 

-1-50° 32' 

’ 44 " 

1.1 

9.1 

0.2800 

3.6 

11.68 

N 

2 

62.7 

ai-slfe" 

21.4+1“/ 

— 0.5 

00 nF27 . 0 
38.0-11 g 

A1763 

13'*35"*18?7 

-1-40° 59' 

' 59 '^ 

0.9 

10.0 

0.2279 

6.9 

14.93 

N 

R 

45.6 

13.6+1;^ 

13.6;“:“ 

36'6-?0.0 

A1835 

14'*01"*02?0 

-1-02° 52' 

'40^^ 

2.2 

14.8 

0.2528 

14.7 

38.53 

Y 

R 

11.5 

o-elo.l 

i-otS'.; 

i2.5;ii 

A1914 

14'*26"*00?8 

-l-37°49' 

'37^^ 

1.0 

10.8 

0.1712 

15.0 

18.39 

N 

2 

29.6 


14 . 2 ;// 

35.9;“2“4 

A2111 

15'*39"*41?1 

-1-34° 25' 

’ 10 ^^ 

2.0 

8.8 

0.2290 

5.0 

10.94 

N 

R 

50.1 

21 . 4 +f^i 

21 . 4 ;“+ 

38.7;«:| 

A2204 

16'*32"*47?0 

-1-05° 34' 

’ 32 ^^ 

5.3 

11.4 

0.1524 

21.9 

21.25 

Y 

R 

7.7 

Q o + O’O 

^•^-0.1 

i7.9;i“ 

A2218 

16'*35"*52?0 

-1-66° 12' 

'34'^ 

3.4 

6.7 

0.1710 

7.5 

9.30 

7 

R 

34.0 

on q+10.7 
'^U.9_i2.o 

30 Q+IO- 7 
'^^■^-12.0 

18 . 4 ;“:/ 

A2219 

16'*40"*20?2 

-|-46°42' 

'30'' 

1.7 

11.4 

0.2281 

9.5 

20.40 

N 

R 

41.1 

29.1+“° 

— 4.5 

29.1+1“ 

— 4.5 

24.9+// 

— 5.6 

A2259 

17'*20"*08?4 

-|-27°40' 

' 10 " 

3.8 

7.1 

0.1640 

5.9 

6.66 

N 

R 

28.9 

^■'^-2.8 

9.o;“:I 

1 c 0 + 8.7 
^^■'^-3.5 

A2261 

17'*22"*27?0 

-1-32° 07' 

' 57 " 

3.2 

10.8 

0.2240 

8.7 

18.18 

N 

R 

18.6 

'^■‘^-0.9 

5 . 0 ;“;“ 

r\cy 0 + 23.5 
^^■'^-8.4 

A2294 

17'*24"*05?5 

-l-85°53' 

’15" 

6.1 

7.1 

0.1780 

5.0 

6.62 

Y 

R 

31.8 

13 1 + 15_5 
— 4.5 

10 1+10.5 
^'^•^-4.5 

64.5;?,“„ 

A2390 

2l'*53"*36?9 

-|-17°41' 

’46" 

6.6 

11.6 

0.2329 

9.6 

21.44 

Y 

R 

20.8 

-^■^-0.3 

1 9+0.4 
^•^-0.3 

ii.7;i:; 

Hercules A 

16'*51"*08?2 

-1-04° 59' 

’ 33 " 

6.0 

6.6 

0.1540 

5.6 

5.62 

7 

R 

17.1 

1 o + O'l 
-^■'^- 0.1 

1 o+O-l 

22.4+1° 

RXJ0439.0+0520 

04'*39"*02?4 

-1-05° 20' 

’ 43 " 

9.7 

7.9 

0.2080 

4.8 

8.62 

Y 

R 

7.8 


4 1 + 2-2 
- 1.1 

19.0+11 

RXJ0439.0+0715 

04'*39"*00?5 

-1-07° 16' 

' 00 " 

1.0 

9.5 

0.2443 

5.3 

13.25 

N 

R 

30.6 

6 1+^'® 

6 1+^-^ 

^ - 1.0 

20,0+9“ 

RXJ1532.9+3021 

15'*32"*53?8 

-1-30° 20' 

’ 59 " 

2.1 

12.5 

0.3615! 

5.0 

24.40 

Y 

R 

10.1 


- 0.1 

15 . 2 ;“'“ 

RXJ1720.1+2638 

17 '* 20"* 10? 1 

-1-26° 37' 

29 " 

3.9 

10.2 

0.1640 

14.3 

16.12 

Y 

R 

13.6 

-^■^-0.3 

1 Q + 0.3 
^•^-0.2 

1 K 9’l"^'® 
^^■^-1.9 

RXJ2129.6+0005 

2l'*29"*39?9 

-1-00° 05' 

’18" 

4.3 

11.0 

0.2350 

8.1 

18.59 

Y 

R 

13.1 

^ - 0.1 

0 9+1 - 0 
'^•^- 0.8 

20 . 7 ;// 

Z1953 

08'*50"*07?0 

-1-36° 04' 

’ 20 " 

3.1 

14.5 

0.3737 

6.1 

34.12 

7 

R 

33.3 

28.8j:+“4 

28.8;^/4 

23.4+1°^’^ 

Z2701 

09'*52"*49?2 

-l-51°53' 

’05" 

0.9 

8.7 

0.2140 

5.6 

10.68 

Y 

R 

13.4 

1 rt + 0.3 
-*-■^-0.1 

2.9j:“-“ 

19.2tt,l 

Z3146 

10'*23"*39?6 

-1-04° 11 ' 

’ 12 " 

2.7 

12.8 

0.2906 

7.7 

26.47 

Y 

R 

11.4 

0.6;“:“ 


12.0t\° 

Z5247 

12'*34"*21?8 

-|-09°47' 

03" 

1.7 

8.5 

0.2290 

4.6 

10.12 

N 

2 

95.5 

169.0;fi,“4 


loo.otfC^ 

Z7160 

14'*57"*15?0 

-|-22°20' 

’34" 

3.1 

9.6 

0.2578 

4.8 

13.19 

Y 

R 

9.3 

0 . 7 ;“:“ 

l-SlS'.a 

21 . 2 ;// 


Col. 1: Cluster name. *The ROSAT derived temperature, flux and luminosity given for A0115 were measured for the source as a whole, rather than as two 
separate, merging clusters; given their projected separation, we treat them as individual clusters in our analysis. Col. 2: Cluster J2000 centroid position 
as determined from the Chandra image after excising contaminating point sources. Col. 3: Galactic absorption column, in units of 10^^ cm~^. Col. 4: 
BCS temperature as derived from Ly^—T relation, in units of keV. Col. 5: BCS redshift. Col. 6: BCS unabsorbed X-ray flux in the 0.1—2.4 keV band, in 
units of 10~ erg~ s~ c m~ Col. 7; BCS in trinsic X-ray luminosity in the 0.1-2.4 keV band (cluster rest frame), in units of 10"^ erg” s. Values in 
columns 3—7 are taken fromexcept for the redshift of RX J1532.9+3021 denoted by “!”, which was determined more accurately by 
Col. 8: Detect ion of Ha line emission as determined byY=yes, N=no, ?=not observed. Col. 9: Cluster 
morphology. Assessed by eye. R=regular, D = Disturbed, 2 = Double peaked, often disturbed. Col. 10: Central cooling bin radius at which the central cooling 
time is measured, in units of kpc. Col. 11: Central cooling time, in units of Gyr. Col. 12: Cooling time at 50 kpc, in units of Gyr. Taken from the closest 
bin below 50 kpc, or interpolated if the central bin is above 50 kpc. (see § details). Col. 13: Cooling time at 180 kpc, in units of Gyr. Taken from 

the closest bin below 180 kpc (see § for details). 


average background rates consistent to better than 15 per 
cent with ACTS quiescent-background calibration measure¬ 
ments.^ 

For such observations, the blank sky backgrounds of 
M. Markevitch, processed in an identical manner to the 
original sample datasets, were used to characterise the X- 
ray background.® For the remaining observations with high 
background rates (noted in column 7 of Table^by “!” or ; 
the symbols respectively denote observations used for both 


^ See http://cxc.harvard.edu/contrib/maxim/bg/index.html 
® See http://hea-www.harvard.edu/~maxim/axaf/acisbg/ 


imaging and spectral analysis and imaging analysis only), 
local backgrounds were determined from regions as free of 
cluster emission as possible. The sample clusters are all small 
enough that we easily can extract source free background 
regions on the 10-13 chips for ACIS-I. However, we caution 
that background regions on the S3 chip for ACIS-S may still 
include some faint cluster emission, thereby affecting the 
spectral analysis for a few of the outermost regions/annuli. 
Note that the degree of such background contamination is 
thus too small to significanly bias any of our conclusions. 

Observations were combined for imaging and source 
detection purposes, using the 0.3-7.0 keV band. Contam¬ 
inating point sources were found using WAVDETECT and 
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Table 1. Summary of Chandra Observations for the z > 0.15 BCS 
Sample. 


(1) 

Cluster 

(2) 

ObsID 

(3) 

Date 

(4) 

ACIS 

(5) 

Exp. 

(6) 

Bgd 

A68 

3250 

2002-09-07 

I 

9.9 

1.039 

A115a 

3233 

2002-10-07 

I 

50.0 

1.045 

A115b 

3233 

2002-10-07 

I 

50.0 

1.045 

A267 

1448 

1999-10-16 

I 

7.8 

1.052 

A520 

528 

2000-10-10 

I 

9.3 

1.010 

A586 

530 

2000-09-05 

I 

9.8 

1.148 

A665 

3586 

2002-12-28 

I 

30.1 

1.068 


531 

1999-12-29 

I 

8.8 

1.126 

A697 

4217 

2002-12-15 

I 

19.7 

0.955 


532 

1999-10-21 

I 

4.4 

1.765* 

A750 

924 

2000-10-02 

I 

29.1 

1.123 

A773 

533 

2000-09-05 

I 

11.1 

1.071 

A781 

534 

2000-10-03 

I 

9.9 

1.032 

A963 

903 

2000-10-11 

S 

36.3 

0.991 

A1204 

2205 

2001-06-01 

I 

23.9 

0.984 

A1423 

538 

2000-07-07 

I 

9.5 

1.019 

A1682 

3244 

2002-10-19 

I 

9.6 

1.784! 

A1758a 

2213 

2001-08-28 

S 

58.6 

1.516! 

A1763 

3591 

2003-08-28 

I 

19.2 

0.935 

A1835 

495 

1999-12-11 

s 

19.7 

1.327! 


496 

2000-04-29 

s 

10.6 

1.323* 

A1914 

542 

1999-11-21 

I 

7.8 

1.329! 

A2111 

544 

2000-03-22 

I 

10.1 

1.084 

A2204 

499 

2000-07-29 

s 

10.1 

1.679! 

A2218 

1666 

2001-08-30 

I 

46.7 

1.010 


1454 

1999-10-19 

I 

11.4 

1.245* 


553 

1999-10-19 

I 

5.7 

1.219* 

A2219 

896 

2000-03-31 

s 

42.5 

1.573! 

A2259 

3245 

2002-09-16 

I 

9.9 

0.990 

A2261 

550 

1999-12-11 

I 

9.1 

1.310! 

A2294 

3246 

2001-12-24 

I 

9.8 

1.106 

A2390 

500 

2000-10-08 

s 

9.6 

1.038 


501 

1999-11-05 

s 

8.8 

1.344* 

Hercules A 

1625 

2001-07-25 

s 

14.8 

1.167! 

RXJ0439.0-1-0520 

527 

2000-08-29 

I 

9.6 

1.003 

RXJ0439.0-1-0715 

3583 

2003-01-04 

I 

18.7 

0.948 


1449 

1999-10-16 

I 

6.0 

1.132 

RXJ1532.9-1-3021 

1649 

2001-08-26 

s 

9.3 

1.056 


1665 

2001-09-06 

I 

9.9 

1.020 

RXJ1720.1-1-2638 

4361 

2002-08-19 

I 

25.7 

0.990 


3224 

2002-10-03 

I 

23.8 

1.035 


1453 

1999-10-19 

I 

7.8 

1.291* 

RXJ2129.6-1-0005 

552 

2000-10-21 

I 

9.9 

0.942 

Z1953 

1659 

2000-10-22 

I 

24.9 

0.971 

Z2701 

3195 

2001-11-04 

s 

27.0 

1.563! 

Z3146 

909 

2000-05-10 

I 

46.4 

0.987 

Z5247 

539 

2000-03-23 

I 

9.1 

1.116 

Z7160 

543 

2000-05-19 

I 

9.6 

1.055 


Col. 1: Cluster name. Col. 2: Chandra observation number. All obser¬ 
vations are used for image analysis, while only the first observation is 
used for spectral analysis. Col. 3: Observation date. Col. 4: Whether 


masked out of the event lists using their 95 per cent 
encir cled-energy radii (e.g., Feigelson, Broos, & Gaffney 
2nnn: rierins et a.i1l2nnnl : M. Karovska and P. Zhao 2001, pri¬ 
vate communication).^ Cluster centres were determined us¬ 
ing a simple centroiding method on the masked 0.3-7.0 keV 
X-ray images (these are given in Table |2J . Radial proHles 
for each cluster were ht initially with single 1-D beta mod¬ 
els, followed by double 1-D beta models if the fit resulted 
in a > 2.0. The core radius, beta index, and amplitude 
were allowed to vary over sensible limits. Postage stamp 0.3- 
7.0 keV X-ray images and radial proHles for each cluster are 
provided in Figure 13 

From inspection of the postage stamp images and ra¬ 
dial profiles in Figure |3 the majority of the clusters appear 
relatively regular and relaxed, although a notable minor¬ 
ity exhibit clear asymmetric morphologies suggestive of re¬ 
cent merger activity or sharp discontinuities indicative of 
inhomogeneities (e.g., cold fronts or shocks). We calculated 
power ratios for each cluster to quantify obvious morpholog¬ 
ical structure within the sample in an objective manner. 

3.2 Power Ratios 

Power ratios allow comparisons of cluster morphologies 
on a common scale, effectively normalising cluster fluxes 
within a spec i fied a perture. Following the prescription in 
iBuote fc Tsail we determined power ratios for each 

of the 38 clusters in the sample. As described therein, the 
power ratios are derived from the multipole expansion of the 
two-dimensional gravitational potential owing to matter in¬ 
terior to a user-defined aperture, essentially measuring the 
square of the ratio of higher order multipole moments to 
the monopole moment (i.e., Pi/Pq, P 2 /P 0 , Ps/Po---)- The 
coordinate system is chosen such that Pi is trivially zero, 
and thus the power ratios P 2 /P 0 , P 3 /P 0 , and P 4 /P 0 provide 
the most useful informatio n for elucidating bas ic morpho¬ 
logical/evolutionary trends iBuote fc Tsailll99fil . e.g.,). Low 
values of Pm/Po indicate highly relaxed, compact and sym¬ 
metric clusters on the scale of the aperture used, while high 
values indicate either relatively flat, spread out clusters or 
clusters with significant substructure or disturbed morpholo¬ 
gies. Odd multipole ratios yield the largest differences be¬ 
tween single and bimodal clusters, essentially vanishing for 
single-component clusters. P 3 /P 0 is particularly sensitive to 
unequal-sized bimodals, with higher values for stronger bi¬ 
modal asymmetry. Even multipole ratios are also able to sep¬ 
arate bimodals from low-ellipticity single-component clus¬ 
ters, but demonstrate more considerable overlap between 
flatter single-component clusters and moderately bimodal 
clusters than the odd terms. 

The power ratios for the 38 BCS clusters are plotted 
in Figure |3 as well as listed in Table |3 We used intrin¬ 
sic physical apertures of 0.5 Mpc to allow comparisons with 


observation was ACIS-I or ACIS-S (i.e., cluster centre was positioned 
on chip 13 or S3, respectively). Col. 5: Useful exposure time after the rnore nearby cluster samples. Images were normalized by 
data were cleaned and background flares removed, in units of ks. Col. 0.3—7.0 keV exposure maps and the average background val- 
6: Ratio of the average 0.3-10 keV background count rate for each ues listed in Tabled were subtracted off. Errors were deter- 
observation versus the ACIS quiescent-background calibration values mined through 1000 Monte Carlo simulations whereby we 
(i.e., 0.31 cte s ^ per ACIS-I chip and 0.86 cts s ^ per ACIS-S3 chip; added Poisson noise to each cluster image and recalculated 
see Sect. 1^ . Symbols indicate that the observation has a high back¬ 
ground, but was used only for image analysis (“*”) or for both image 

and spectral analysis ( ! ). 4 Peigelson et al. (2000) is available at 

http: / / www.astro.psu.edu / xray / acis / memos / memoindex.html 
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Figure 2. Almages are too large to include in astro-ph version; download full version at 

http://www.astro.Columbia.edu'^feb/Bauer.ps.gz Left: S' x 8^ colour images for each BCS cluster in the sample with Chandra 
data, comprised from the 0.5-1.0 keV (red), 1-2 keV (green), and 2-7 keV (blue) X-ray energy bands. The images have been spatially 
binned by a factor of 2, smoothed with a a =3 pixel Gaussian smoothing scale, and logarithmically scaled (following the methods of 
iLupton et alJl2004f) . The white cross denotes the cluster centroid, while the white circle indicates a radius of 0.5 Mpc at the distance of 
each cluster. Middle: Radial profiles for each BCS cluster in the sample. Right: Projected temperature maps for each BCS cluster in the 
sample, scaled between 0—20 keV. The maps have been cropped to the 0.5 Mpc radius shown in the colour images to speed computation. 
The signal-to-noise of the data in each region (which varies between 20—40) is shown in the bottom right corner of each map. 


power ratios. As zero is a very poor estimate of the Poisson 
noise, we instead adopted the average background value per 
pixel (typically quite small) for each empty pixel of a given 
observation. 

We Hnd that the power ratios from this sample oc¬ 
cupy nearly identical paramete r space to those fro m the 
nearby sample of clusters from iBuote fc Tsail (Il99(jl . The 
only deviations are from the extremely compact clusters 
A0115a, A2294, and RXJ0439.0-1-0520 which extend to very 
low PijPg values. Overall, the majority of clusters in our 
sample he in the region of power ratio parameter space typ¬ 
ical of relaxed, symmetric clusters. 


3.3 Spectral Analysis 

Spectral analysis of the clusters was performed following 
two methods. For simplicity. X-ray spectral analysis in both 
cases was always performed only on the hrst entry in Tabled 
(typically the longest observation) and data were masked out 
beyond intrinsic radii of 0.5 Mpc to simplify and speed up 
analyses. 

In the first method, we divided the da ta into several 
regio ns using a “contour binning” method llSanders et alJ 
I 2 OO 5 I : 2005b, in preparation) such that the signal-to-noise 
was 20-40 per region depending on the quality of the data. 
The algorithm defines regions with a signal-to-noise greater 
than a threshold by growing bins in the direction on a 
smoothed map which has a value closest to the mean value 
of those pixels already binned. This technique dehnes bins 
which are matched to the surface brightness distribution of 
the object. X-ray spectra were extracted for each region and 
fit in XSPEC with a single-temperature MEKAL model ab¬ 
sorbed by a PHABS model. The temperature in each region 
was allowed to vary between 0.1-20 keV, while the absorp¬ 
tion was fix ed at the Galactic v alue along the line of sight 
(taken from lEbeling et al.lll99s|) . the abundance was fixed 
to 30 per cent of solar,® and the redshift was set to its ap¬ 
propriate value for each cluster. The results of this spectral 
analysis were combined together to form the projected tem¬ 
perature maps shown in Figure]^ The projected tempera¬ 
ture maps complement the colour images in terms of iden¬ 
tifying clusters with strong cool cores, but are particularly 
useful for revealing weak cool cores and further substructure 
not readily apparent from the images alone. 

In the second method, we split the data into no more 
than 10 concentric circular annuli again such that the signal- 
to-noise was at least 20 per region, but often much higher 
depending on the quality of the data.® X-ray spectra were ex- 

® Letting the abundance vary makes little difference to the mea¬ 
sured temperatures. 

® Most of the clusters in our sample are in fact elliptical rather 


tracted for each region and these regions were fit simultane¬ 
ously in XSPEC with a PROJCT model, adopting absorbed 
(PHABS) single-temperature MEKAL models for each an¬ 
nular bin (i.e., a model identical to that used in the first 
method). Using these annular spectral fits, we derived depro- 
jected radial temperature, density, cooling time, and mass 
deposition rate profiles. 


4 DISCUSSION 

As found in Sections 13. ll and |3.21 the majority of the clusters 
have relatively regular morphologies and generally appear to 
be relaxed. From a visual inspection of the 38 BCS clusters, 
we find five that are obviously double-peaked and four more 
that have disturbed morphologies indicative of merging (see 
Table H The temperature maps in a handful of clusters hint 
at mild additional substructure within 0.5 Mpc. We note 
that our assumption of spherical symmetry does little to ac¬ 
count for clusters with double-peaked or merging morpholo¬ 
gies, which appear to make up « 25 per cent of the sample. 
In these cases, our radial profiles simply serve as a proxy 
and are likely to overestimate the real temperature and un¬ 
derestimate the real central cooling time. Thus we caveat 
that cooling times for disturbed clusters should be regarded 
as likely lower limits. We still provide this analysis since 
it may be useful to have some point of reference for com¬ 
parison to lower-resolution nearby clusters observed with, 
e.g., ROSAT or high-redshift clusters observed with XMM- 
Newton or Chandra, where statistics are generally poorer 
and mergers less obvious. 


4.1 Cooling Time Profiles 

We have measured the cooling time profiles for all of the 
clusters in our sample, as shown in Figure |1] While the 
cooling times in the central bin differ by a few orders in 
magnitude, we note that there appears to be a prominent, 
apparently universal lower bound to all of the cooling pro¬ 
files. The average cooling time profile for the 13 clusters with 
a central cooling time less than 2 Gyr (thick black line) in¬ 
dicates that there is little difference in the lower bounds for 
a set of clusters with widely ranging physical characteristics 

than spherical. However, for simplicity we have derived profiles 
using circular rather than elliptical annuli since we do not know 
whether the clusters are prolate or oblate. To determine the de¬ 
gree of error this might introduce, we fitted elliptical annuli to a 
few of the most eccentric clusters in our sample. We found virtu¬ 
ally no change in the derived profiles; the circular and elliptical 
curves were consistent with each other to at least a ~ Icr confi¬ 
dence level with no apparent systematics. 
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Figure 3. Power ratios for the 38 BCS clusters in the sample using a 0.5 Mpc aperture. Filled black, grey, and empty symbols denote 
clusters with strong cool cores {tc < 2 Gyr), weak cool cores {tc = 2-10 Gyr), and no obvious cool cores (tc > 10 Gyr), respectively. 
Gircles denote morphologically regular clusters, while triangles denote double-peaked or disturbed clusters. A2294 (open circle in the 
lower left of the first two panels) is again an exception, being a very compact cluster without a cool core. 



Figure 4. Radial cooling time profiles for each of the 38 BGS 
clusters in the sample (denoted by the different grey/colour lines, 
for paper/online publishing respectively). There appears to be a 
lower bound to the cooling profiles. The thick black line, which 
traces this lower bound, denotes the average cooling profile for 
the clusters with a central cooling time less than 2 Gyr. 

(i.e., temperatures, densities, entropies). While we might ex¬ 
pect some sort of boundary given the general nature of how 
gas cools in clusters, it is physically unclear why a univer¬ 
sal boundary that takes the form tc oc (depending 

on w hether one imposes an inner cutoff) would exist. Simi¬ 
larly, IVoigt fc FabiarJ ^2^^4^ find a very small dispersion in 
cooling time profiles among their sample of 16 clusters with 
strong cool cores, with tc oc . 

The cooling times for the innermost radial bin are shown 
in Figurel^for each cluster. Unfortunately in physical units, 
the innermost radial bin is both a function of the cluster 
distance and the overall signal-to-noise of the observation 
(which in turn depends on whether a cluster exhibits a cool 
core); hence the apparent correlation in Figure]^ Impor¬ 
tantly, the spatial resolution of Chandra is a factor of several 
smaller than the typical innermost radius adopted (i.e., 0/5 
corresponds to 1-3 kpc), so any failure to detect a significant 
cool core in the BCS sample is not due to any limitation on 


Figure 5. Plot of the central cooling time for the inner radial bin 
for the 38 BCS clusters in the sample. Filled black, empty, and 
filled grey symbols denote clusters with det ected, undetected, an d 
unobserved Ha line emission, respectively j£rawfordjet_ak|l22SD • 
Circles denote morphologically regular clusters, while triangles 
denote double-peaked or disturbed clusters. There is a strong cor¬ 
respondence between low cooling time and detection of Ha line 
emission. The one exception is A2294. 

spatial resolution. Also note that because we use the cluster 
centroid rather than the surface brightness peak, the cool¬ 
ing times inside a few tens of kpc could be biased upward by 
at most a factor of Alternatively, the central cooling 

times should be considered upper limits because the strength 
of cooling depends strongly on the radius at which it is mea¬ 
sured, and thus on the signal-to-noise of the observation. 

From our determination of cooling times for the inner- 


^ The peak and centroid surface brightnesses (oc F) and clus¬ 
ter temperatures never differ by more than a factor of 2 each. 
Since tc oc T® ® n ^ , this means that tc is unlikely to change dra¬ 
matically based on our adopted position. Significant differences 
between the peak and centroid positions only occur in highly dis¬ 
turbed or double-peaked clusters and comprise at worst < 20 per 
cent of our sample (see Table O- 
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Table 3. Power Ratios for the Cool Core Sample. 


(1) 

(2) 



(3) 



(4) 

Cluster 

P 2 /P 0 


P 3 /P 0 


P 4 /P 0 

A68 

(i-2itrii) 

X 

10-® 

(1.02l°;«°) 

X 

10-7 

( 5 . 8 OII!!) X 10“® 

A115a 

(2-40i;,i^) 

X 

10-7 

(1.86in?) 

X 

10-7 

(1.1811??) X 10-7 

A115b 

(3.72+°;“) 

X 

10-® 

(3.871°;°!) 

X 

10-7 

(1.0911?!) X 10-7 

A267 

(7.23li;ll) 

X 

10-° 

(3.1411™) 

X 

10-® 

(2.251?-™) X 10-8 

A520 

(6.08li;“) 

X 

10-° 

(2.88llt^) 

X 

10-7 

(1.7611™) X 10-7 

A586 

(8.18+2-2®) 

X 

10-7 

(2.4511™) 

X 

10-® 

(1.561!-™) X 10-8 

A665 

(3.27t°-®2) 

X 

10-° 

(1.64ll^^) 

X 

10-7 

(1.3311?!) X 10-7 

A697 

(6.55t°;S) 

X 

10-° 

(1.56ll«t) 

X 

10-7 

(1.471!;!?) X 10-« 

A750 

(2.33+°-28) 

X 

10 -° 

(4.25ll°l) 

X 

10-8 

(7.6811™) X 10-!° 

A773 


X 

10 -° 

(2.9411!°) 

X 

10-8 

( 6 . 97 I?-??) X 10-8 

A781 

(4.67t+,f5) 

X 

10 -° 

(1.44ll^?) 

X 

10-7 

(4.181?-??) X 10-7 

A963 

(7-00l+,“) 

X 

10-7 

(s.iilli^) 

X 

10-8 

(1.511!-??) X 10-° 

A1204 

(3.72«;®®) 

X 

10-7 

(1.1811!!) 

X 

10-8 

( 6 . 91 I?;™) X 10-° 

A1423 

(3.65«-®4) 

X 

10 -° 

(3.5911!!) 

X 

10-8 

(1.961? !?) X 10-® 

A1682 


X 

10 -® 

(1.3611!!) 

X 

10 -° 

( 3 . 37 I?-?!) X 10-7 

A1758a 

(4.20l°;«) 

X 

10 -® 

(2.4611™) 

X 

10-7 

(4.6611™) X 10-7 

A1763 

(9.74t°-®®) 

X 

10 -° 

(9.631!;!!) 

X 

10-8 

(3.671?-?!) X 10-8 

A1835 

(6.26+°;®!) 

X 

10-7 

(2.3811™) 

X 

10-8 

( 6 . 72 I!;™) X 10-™ 

A1914 

(1.94l°;24) 

X 

10 -° 

(1.1711!?) 

X 

10-7 

(1.12l!-7°) X 10-!° 

A2111 

(i.45t°;?®) 

X 

10 -® 

(1.4711?!) 

X 

10-7 

(1.041?-™) X 10-7 

A2204 

(5.32+°;«) 

X 

10-7 

(1.0511!!) 

X 

10-8 

(7.641?-??) X 10-° 

A2218 

(2.65l°-2°) 

X 

10 -° 

(i.isllw) 

X 

10-7 

(1.981?-™) X 10-° 

A2219 

(2.081°;°°) 

X 

10 -® 

(I. 61 II!!) 

X 

10-7 

(1.3911?!) X 10-7 

A2259 

(4.081°;™) 

X 

10 -° 

(1.2811™) 

X 

10-7 

(1.041!-™) X 10-8 

A2261 

(7.4211^2) 

X 

10-7 

(4.1411™) 

X 

10-8 

(1.191!;?!) X io-« 

A2294 

(4.9611,!^) 

X 

10 -® 

(6.7611™) 

X 

10 -° 

( 8 .I 2 I!-??) X 10-1° 

A2390 

(7.12l°-J°) 

X 

io-« 

(2.1511!?) 

X 

10-8 

(3.091?-??) X 10-8 

Hercules A 

(i.i8l°;“) 

X 

10 -° 

(6.831!-!!) 

X 

10 -° 

(4.911?;!?) X 10-1° 

RXJ0439.0-I-0520 

(8.6lll,°l) 

X 

10 -® 

(7.1311?!) 

X 

10 -° 

(1.131! !?) X 10-° 

RXJ0439.0-1-0715 

(3.32l°-«) 

X 

10 -° 

(i.isllS) 

X 

10-7 

(1.251?-??) X 10-8 

RXJ1532.9-1-3021 

(6.0711™) 

X 

10-7 

(1.061!;!!) 

X 

10 -° 

(4.281?-??) X 10-1° 

RXJ1720.1-1-2638 

(9.77ll“) 

X 

10-7 

(5.8911??) 

X 

10 -° 

(4.571?-??) X 10-° 

RXJ2129.6-I-0005 

(1.871°;°°) 

X 

10 -° 

( 4 . 61 II™) 

X 

10 -° 

(2.661?-??) X 10-8 

Z1953 

(4.021°;™) 

X 

10 -° 

(1.6411!!) 

X 

10-7 

(I. 77 I?-??) X 10-8 

Z2701 

(i.07l°;°°) 

X 

10 -° 

(8.5i1!°5!) 

X 10-!° 

(4.131?-!?) X 10-1° 

Z3146 

(6.68l°;g) 

X 

10-7 

(6.6811™) X 10“° 

(2.051!;!?) X 10-° 

Z5247 

(7.0il°;f2) 

X 

10 -® 

(1.941?-™) X 10-7 

(1.491?-!?) X 10-° 

Z7160 

(3.9911°°) 

X 

10-7 

(3.8911™) X 10-° 

(3.331?-?!) X 10-° 


Col. 1: Cluster name. Alternate names are given in brackets. Cols. 2—4: Power ratios 
as determined in Sect. 13.21 for an aperture of 0.5 Mpc. 


most radial bin, we find that at least 20 of 38 (53 per cent) 
BCS clusters in the Chandra sample have central cooling 
times < 10 Gyr (mild cool cores) and at least 13 of 38 (34 
per cent) have < 2 Gyr (strong cool cores). To test for any 
evolution, we compare the cumulative cooling time fractions 
from the BCS sample to those found for the nearby B55 sam¬ 
ple of IPeres et alJ (1199811 , who provide cooling times for both 
the innermost radial bin and 180 kpc (originally 250 kpc in 
their cosmology). At first glance, the fraction of cool core 
clusters in our sample as measured from the innermost ra¬ 
dial bin would appear to be markedly different from those 
in the B55 sample (see left plot in Figure]^, implying sub¬ 
stantial evolution in the cores of X-ray bright clusters out 
to 2 ~ 0.4. However, we must exercise caution here, since a 
direct comparison between the two samples in this manner 


is problematic. Importantly, the innermost radial bins span 
a large range in both samples and are significantly smaller 
in the B55; thus comparing these quantities is not likely 
to provid e a useful yardstic k for cooling. Additionally, we 
note that IPeres et alJ ^1998^ measured their cooling times 
using the surface brightness peak rather than the centroid, 
so there may be a slight systematic difference between the 
calculated cooling times for a small fraction of clusters. A 
final consideration is that the ROSAT observations of the 
fainter or more distant clusters in the B55 sample may fail 
to adequately resolve potential cool cores due to limited spa¬ 
tial resolution. 

Given the above problems, we instead choose to com¬ 
pare cooling times for the two samples at a fixed intrinsic 
radius where cooling is still important. Since a large frac- 
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tion of both samples have innermost radial bins close to or 
smaller than 50 kpc, this value seems like a sensible radius 
to use. At 50 kpc, we find that at least 21 of 38 (55 per 
cent) BCS clusters in the Chandra sample have central cool¬ 
ing times < 10 Gyr (mild cool cores) and at least 8 of 38 (21 
per cent) have < 2 Gyr (strong cool cores). Because the B55 
sample lacks published cooling time profiles, we must inter¬ 
polate cooling times at 50 kpc. For this purpose, we adopt as 
an interpolation template the average cooling profile found 
in Figure |1] for clusters with tc < 2 Gyr. Comparisons at 
50 kpc (middle), as well as at 180 kpc (right), are shown 
in Figure El The apparent differences in the cooling time 
fractions measured using the innermost radial bin have now 
all but vanished, demonstrating that there is no detectable 
evolution in the cool core fraction of clusters out to ^ ~ 0.4. 
Note that the small discrepancies for short cooling times 
may be due to our assumed interpolation template. Also 
note that the fraction of cool cores in either sample has no 
obvious X-ray luminosity dependence. 

Let us consider the evolution of a cluster core host¬ 
ing a cooling flow. First assume that it is isolated with no 
heating, so that radiative cooling proceeds to drop the gas 
temperature. Approximating the situation locally as a con¬ 
stant pressure flow, then the cooling time at 50 kpc will 
decrease by approximately the mean time difference be¬ 
tween the samples. The mean redshift of the B55 sample 
is 0.056 and of the 38 BCS clusters is 0.22, so this differ¬ 
ence in our adopted cosmology is about 1.8 Gyr. It is clear 
from Fig. El that this is not occurring since at the 50 per 
cent level the curves are consistent with each other. Second, 
we consider the effect o f the continual accretion of sub clus- 
ters by a cluster (e.g., iR.owlev. Thomas, fc Kavll2004^ . so 
that the cool central region undergo adiabatic compression 
which, for bremsstrahlung cooling appropriate here, varies 
as T~^. This again reduces as time proceeds. Weak to mod¬ 
erate shocks also reduce the cooling time. What is needed 
is a non-gravitational heat source which offsets the effects 
of radiative cooling. The similarity between the cumulative 
curves at 50 kpc for the B55 and BGS samples demonstrates 
that the heating invoked to balance radiative cooling re¬ 
quired to explain the detailed temperature profiles of the 
clusters must have been in place before the redshift of the 
BCS. More distant samples are required to determine just 
when the heating/cooling balance was established. 

4.2 Ubiquity of a Cooling “Floor” 

As noted in § Q the central temperatures in clusters are 
found to generally drop by less t han a factor of ~3- 
4 with very l ittle c ooler gas (e.g., |Petepoi^^lJj20Ml 
[r ^jmira et al.l I 2 OOII : IPeterson et al'l 120031 : iKaasto^e^^ 
1222^. This cooling “floor” is thought to be due to some 
form of heating (e.g., from an active galactic nucleus, con¬ 
duction). Figure El compares the temperature drop — taken 
to be ratio between the central bin temperature and the av¬ 
erage temperature outside of 180 kpc — to both the 0.5 Mpc 
power ratios and central cooling times (right) for the 38 BCS 
clusters in the sample. A similar cooling floor (i.e., a tem¬ 
perature drop of no more than « 3-4) exists for our cluster 
sample, and does not appear to correlate with either the 
strength of the cool core or its general morphology. Inter¬ 
estingly, three of the four lowest temperature drop clusters 


are merging clusters; inspection of the temperature maps for 
these clusters indicates that their large temperature ratios 
are in part likely due to shocks in the outer regions of those 
clusters. In general, we find that the disturbed clusters in 
our sample, which all have high P3/P0 values, span the en¬ 
tire range of observed temperature ratios, presumably due 
to the varied merger histories and ages of these clusters. 


4.3 Cooling Time vs. Ho 

A significant fraction of clusters with cool cores have also 
been shown to exhibit signs of filamentary Ha line emis¬ 
sion. Since ne arly all of the clusters in our sample have Ha 
observations iCrawford et alJll99flll . we can examine how 
useful Ha emission is as an indicator of cool cores and clus¬ 
ter morphological properties. FigureEldemonstrates that 13 
of the 38 clusters in the sample show evidence for Ha line 
emission and all but one of these are among the clusters 
with the very shortest central cooli ng times.^ This re sult 
strengthens the correlation found bv IPeres et alJ ^1998^ for 
relatively nearby clusters. Interestingly, there are another 
three clusters among the 18 not observed by Chandra that 
also show evidence for Ha line emission. If the Ha-detected 
fraction traces the overall cluster population, this would sug¬ 
gest that ~15-20 per cent of the unobserved sources are 
likely to have strong cool cores and that we might expect a 
similar overall cool core distribution (i.e., both strong and 
moderate) among the unobserved clusters as found for the 
sample sources. 

The observed Ha and star formation rates requires the 
mass of Ha-emitting gas typically to be only a few to a 
few tens of Mq yr“^, whereas classical cooling flows in the 
observed clusters would require many hundreds of Mq yr“^. 
Thus the presence of Ha could imply continuous cooling 
below the “floor” at a rate of up to a few tens of Mq yr“^. 
Or alternatively, it could imply the presence of conditions 
necessary for long-lived Ha. 


4.4 Comparison to Power Ratios 

A large number of higher redshift clusters have now been 
found. As these sources are less likely to have the signal-to- 
noise or spatial resolution to unequivocally detect any po¬ 
tential cool cores, it is instructive to compare our cool core 
estimates and Ha detections for the BGS clusters with their 
power ratios (which can be measured more reliably under 
the above conditions). Indeed, FigureEldemonstrates that 
the power ratios provide a quantitive assessment of basic 
morphological trends and can be used as a proxy for cen¬ 
tral cooling times. Figure El bears out this relationship bet¬ 
ter, showing the strong est correlation betwe en P3/P0 versus 
tcooi. As expected ('e.g.. lBuote fc Ts^ll99(tl . the more sym¬ 
metric, compact clusters generally have stronger cool cores 
and a higher likelihood of Ha line detections. 


® A2294 is the one outlier. This source has a strong detection 
(C. Crawford, private communication), although we are unable 
to rule out the possibility of contamination from nearby galaxies. 

















The Prevalence of Cooling Cores in Clusters of Calaxies at z ^ 0.15-0.4 9 





Figure 6. Comparison of the cumulative fraction of cooling time s ( left : cent ral bin; middle: at 50 kpc; right: at 180 kpc) for the 38 BCS 
clusters in the sample (black) and the B55 sample of IPeres et al.l fgrev: [l998ll . The dash-dotted lines denote the 1 a error bars. The large 
apparent difference seen in the cumulative fraction measured for the central bin disappears when comparisons are made at fixed intrinsic 
radii. 
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Figure 7. Plot of the ratio between the central bin temperature and the average temperature outside of 180 kpc versus 0.5 Mpc power 
ratios (left) and central cooling times (right) for the 38 BCS clusters in the sample. Left: Filled black, grey, and empty symbols denote 
clusters with strong cool cores (t^ < 2 Gyr), weak cool cores (tc = 2-10 Gyr), and no obvious cool cores (i^ > 10 Gyr), respectively. Right: 
Filled_blackj_empt 2 jand filled grey symbols denote clusters with detected, undetected, and unobserved Ha line emission, respectively 
JCrawford et alJ ll99fl^. Circles denote morphologically regular clusters, while triangles denote double-peaked or disturbed clusters. 


5 CONCLUSIONS 

Using a sample of 38 luminous BCS clusters at « = 0.15-0.4, 
we find that cool cores still appear to be common at mod¬ 
erate redshift. At 50 kpc, at least 55 per cent of the clusters 
in our sample possess mild cool cores (tcooi < 10 Gyr) and 
within the central bin at least 34 per cent possess strong 
cool cores (tcooi < 2 Gyr). We find that a cooling floor ex¬ 
ists for our moderate redshift sample, similar to that found 
for many nearby clusters, such that the ratio of central to 
outer temperatures rarely increases above a factor of ss 3-4. 
Moreover, comparing the central cooling times to catalogues 
of central Ha emission in BCS clusters, we find a strong cor¬ 
respondence between the detection of Ha and short central 
cooling times. We also find a correlation between the cen¬ 
tral cooling time and cluster power ratios, indicating that 
crude morphological measures are a proxy for more rigorous 
analysis in the face of limited signal-to-noise data. 

Comparing our moderate redshift BCS sample to the 
local B55 sample [z < 0.15), we find no evidence for any 
significant evolution in the cumulative cooling time fractions 


at both 50 kpc and 180 kpc. This suggests that cool cores 
are already well-developed even at moderate redshifts and, 
as a consequence, are quite likely to be robust against a 
wide variety of merger scenarios. It also implies that any 
additional cooling our sample of clusters underwent over the 
last ~2-3 Gyr must be balanced by some form of heating. 
Thus, balanced heating and cooling mechanisms are likely 
to have already stabilised or been frozen in by this epoch. 

Looking toward the future, it will be instructive to anal¬ 
yse Chandra and XMM-Newton observations for a sample of 
local clusters (such as the B55 sample) in a manner similar 
to that done here for the BGS sample, so we can perform a 
more complete and rigorous comparison of the two samples. 
Further work is also needed to understand when the onset of 
cool cores occurs. The high fraction of cooling core clusters 
observed here needs to be tied to the relative dearth of such 
clusters at the highest redshifts (z ~ 1). Another issue is to 
determine when and how the heating process, which leads 
to the nearly ubiquitous temperature plateau of a factor of 
3-4, occurs. This heating process appears to already be in 
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Central Cooling Time (Gyr) 

Figure 8. Plot of the P 3 /P 0 power ratios and central cooling 
times for the 38 BCS clusters in the sample. Filled black, empty, 
and filled grey symbols denote clusters with detectedjjmdetecte^ 
and unobserved Ho line emission, respectively JCrawford_e^,ajJ 
|]^99|). Circles denote morphologically regular clusters, while tri¬ 
angles denote double-peaked or disturbed clusters. 

place for the clusters in our sample, so again we must look 
to higher redshift samples for clues. 
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